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Abstract 

The Hirota bilinear difference equation is generalized to discrete space of arbi- 
trary dimension. Solutions to the nonlinear difference equations can be obtained via 
Backlund transformation of the corresponding linear problems. 
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1 Introduction 

Hirota bilinear difference equation (HBDE) |Q plays the central role in the study of integrable nonlinear 
systems. This single equation embodies infinite number of integrable differential equations which belong to 
the KP-hierarchy as shown by Miwaj^j,|Q. It characterizes algebraic curves J3|. It appears as a consistency 
relation for the Laplace maps on the discrete surface H. It is satisfied by string correlation functions in 
particle physics ||. It is also satisfied by transfer matrices of some solvable lattice models || |?J Q |]]. 

Solutions to HBDE, called r function, are known |To| to form a Grassmann manifold of infinite 
dimensions. On this space of solutions a Backhand transformation acts and generates the GL(oo) sym- 
metry. This large symmetry is the origin of integrability of the system. The scheme of the Backlund 
transformation can be seen most easily if we convert HBDE into a pair of linear homogeneous equa- 
tions. It was first found by Hirota in El and then represented in a manifestly symmetric form under 
the Backlund transformations in ref. [fl2| . The coupled linear equations describe the behaviour of a wave 
function under the influence of gauge potentials. The compatibility condition for the wave function to 
solve the coupled equations requires the gauge potential to satisfy HBDE. 

The situation is similar to the inverse scattering method. The difference from the ordinary inverse 
method is that every solution to the linear equations also fulfills the same nonlinear equation satisfied by 
the gauge potential. This fact owes to the remarkable symmetry under the exchange of the wave function 
and the gauge potential, which we called the dual symmetry in ]l2| ]. Hence, by solving the linear system 
of equations iteratively we obtain a series of solutions to the nonlinear equation. We call this scheme of 
finding new solutions to a nonlinear equation, a linear Backlund transformation (LBT). 

An interesting application of this scheme was found and described in The authors showed that 
it can be used to generate Bethe ansatz solutions to certain class of solvable lattice models. In order to 
uncover the meaning of the correlation between the solvable lattice models and the HBDE, it is desireble 
to formulate LBT such that the symmetries of this scheme can be seen from both sides. From this point 
of view we like to discuss, in this paper, a symmetric generalization of LBT. 

The HBDE is a highly symmetric equation under the exchange of lattice variables. The linear equa- 



tions discussed in 12 , however, do not possess this symmetry. In §3 we derive a symmetric version 
of the LBT associated with HBDE. We will show in §4 that we can generalize this scheme to lattice 
space of arbitrary dimension. Corresponding to this generalization we obtain a large number of nonlinear 
equations which can be solved by the LBT method. The HBDE turns out to be a special case of this 
scheme. 
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2 Hirota Bilinear Difference Equation (HBDE) 

Let us begin with writing the Hirota bilinear difference equation (HBDE)|l|] for a function / 6 C°° of 
discrete variables A, /x, v. 

af{\ + l,n,v)f{\-l,ii,v) + /3/(A,m + M/(A,m-1,^) 

+ 7/(A,/x,2. + 1)/(A,m,^-1)=0 (1) 

where a, /3, 7 are arbitrary complex parameters. This is a very simple and highly symmetric equation. 
In order to show how the LBT generates solutions to ([!]), we follow the argument of ref. j5j, but in a 
different notation convenient later. 

We consider the following coupled linear problems: 

Vi 2 9 = W12 5, V 2 i g = OJ21 5, (2) 
V12 := f2e B »- B »fc 1 -CLfie a »- B »fr 1 

V21 := /xe-^-^/r 1 - ca/ae-^-^/a- 1 , (3) 

where /1, /a are functions of A, /x, ^ and W12, W21, ci, C2 are constants. The shift operator e 9a: acts to 
all functions on the right by changing x to x + 1. 

HBDE (|l|) emerges from the compatibility condition of this set of linear equations 

[Vi 2 , V21] = 0. (4) 

In fact, if we impose the condition /i(A, /x + 1, v + 1) = /2(A + l,fi, v + 1) =: /(A, /x, z/) to the gauge 
potentials, this operator relation turns to 

/(A-l,/i-l^)/(A,/i,i/) 16 ' 6 j 
^ f /(A + l, M ,^)/(A-l, M ,^) f(\,y + l,v)f(\,Li-l,v) \ _ 

\f(\,H,v + l)f(\,H,v-l) ClC ' 2 f{\^u + l)f{\^u-l)\ ■ [) 

Since the left hand side is a difference of the quantity in the brace, this quantity itself must be a constant. 
By adjusting the parameters as C1C2 = —/3/a, and putting the constant — 7/a we obtain the HBDE. 

An important observation is that the linear equations (Q) are symmetric under the exchange of the role 
of the gauge potentials / and the wave function g. This owes to the fact that under the same constraints 
for the guage fields /'s, (||) can be also written as 

Via / = wia /, V21 / = wai /, (6) 
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if we define gx(X, n — 1, v — 1) = g 2 {\ — 1, /i, f — 1) = g(A, /i, ^) and 



Via : 
V21 := g\e 



9 ^g^-c 2 g 2 e 9 ^g;\ 



(7) 



It is not difficult to convince ourselves that the linear equations @ for the new wave function / can be 
solved only if the new gauge potential g satisfies (||), the same equation satisfied by /. This is what we 
call dual symmetry in [l^j . 

The linear Backlund transformation works as follows: Starting from a particular solution to HBDE, 
say f( x \ as a gauge potential we solve the linear problem (||). One of its solutions, which we call g^ 
satisfies the HBDE because it can be regarded as a potential of the coupled linear equations @ which 
are nothing but another expression of (j^). 

The set of coupled equations (|6|) should have solutions other than / W. Let us call it f^ 2 \ must 
also satisfy HBDE, because (||) has g^ as a solution, hence it is compatible. Now we can repeat the 
same argument starting from to obtain new solution g^ to HBDE, and so on. This is an auto 
Backlund transformation since /'s and </'s are solutions to the same equation. By starting from the 
simplest solution / W = 1 in the case of a + (3 + 7 = in (|l|), for instance, we obtain a series of soliton 
solutions explicitly in this method. 

3 Symmetrization of LBT 





Fig. 1 



Fig. 2 



The HBDE is highly symmetric by itself, in contrast to its continuous reductions, such as the KdV 
equation, Toda lattice equation, sine-Gordon equation etc.. The corresponding linear version of HBDE, 
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(^) and (Q), however, are not symmetric. In this section, we like to derive complete set of linear equations 
which recover the symmetries possessed by HBDE. 

First we examine the symmetries of the HBDE. From (0) HBDE relates functions defined on the six 
lattice sites (A ± 1, (i, v), (A, fj, ± 1, v), (A, fi, v ± 1), which form an octahedron in the three dimensional 
lattice space. Connecting these corners of octahedrons forms a face centered cubic (fee) lattice. Hence 
dependent variables of HBDE reside on fee lattice, rather than a simple cubic lattice. Indeed the vectors 
p, q, r in Fig.l are the primary translation vectors of fee. 

The HBDE possesses point group symmetries which transform an octahedron to itself. It includes 
the inversion /, which changes the sign of the variables A, /i, v, and the mirror reflections ov,- which 
exchange the triangles Si and Sj in Fig. 2. The linear Backlund transformation represented in the form 
of (^|), however, is not symmetric under the group transformation. The purpose of this section is to make 
it symmetric. 

First we notice that the potential field f\ appears in (|2|) as a gauge field in the form of covariant shift 
operators /le^" 9 " /-f 1 and fie~ dx ~ d " /-f 1 defined on two edges of the triangle S\ in Fig. 2. Similarly the 
field fi is associated with the triangle S2. On the other hand from (^) we learn that the fields g\ and g 2 
are associated with the triangles S\ and S2, respectively, which are obtained from Si and S2 via inversions 
with respect to the origin. Hence two equations (||) and (||) are related each other by the inversion /. 
In other words the inversion is realized by the dual symmetry. The rest of the point group symmetry 
requires potential fields associated with S3 and S4 to exist. Let us call them /a and /4, respectively. 

In order to make the symmetry manifest it is convenient to introduce the following set of derivatives 

d x = , d 2 - , 3 = g ' di 2 (8) 

They represent the gradients along the normal to the surfaces Si, S2, S3, S4. In terms of these operators 
the point symmetry transformations can be expressed as 

on ■ fi^fj> di^dj. (9) 
Applying cr's to (^) successively we obtain, up to coefficients, the following six pairs of equations 
{fy ■' !\ : + fke^fk 1 + l)g = 0, 

(.ny 1 + i/ K "l 1 + 1)9 = 0, 

j < fc, (i,j,k,l) = even permutation of (1,2,3,4). (10) 



Among the twelve equations of (|10|) we notice that, corresponding to each triangle of the octahedron, 
there are three equations which connect two edges of the same triangle by the gauge fields. For instance 
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the set of equations (j,k) = (1,2), (2,4), (1,4) determine the triangle S3. For the consistency of these 
equations they must be the same. We can resolve this problem if we adopt the following conditions to 
the gauge fields: 

/i(A,At,i/)=e*- & /(A,At,i/) ) 3 = 1,2,3,4. (11) 
Using the fact that d\ + 82 + B3 + §4 = we can write ( fici|) as 

(e-h-bf) (e^ +Si g) + (e^*"*/) (e 3 * +8A g\ + (e"^"*/) (e^+S) = (12) 

for all even permutations of (1,2,3,4). Here we used the notation (e 9 /) to mean that e d acts only to the 
functions in the bracket. In this expression it is obvious that, for a given i, three equations with different 
choices of j, k,l are the same equation. Therefore we have only four different equations, which we can 
write as 

5> y - (e-^'/) (<" "-/) =0, i = 1,2,3,4 (13) 

We have recovered coefficients ay in this expression. They are free unless we specify the values of 
a, (3, 7 in HBDE. This is the symmetrized LBT which we were looking for and is symmetric under the 
a transformations. 

The symmetric LBT ( [f3| ) possesses the a symmetry as we required. From the analogy of the cor- 
respondence between (^) and (Q), the inversion symmetry, which has not been imposed so far, should 
be included if the dual symmetry under the exchange of / and g holds. In fact if we multiply the shift 
operator e di ~ di from the left to the ith equation of (|l3|), it becomes 

£ o« (e* + %) (e-^-M =0, i = 1,2,3,4. (14) 

These are exactly the equations we obtain from (|l^) by the inversion transformation / of ([)]). Therefore 
( |l3| ) itself is the totally symmetric LBT which we expected. 

Now we will show how HBDE arises from this set of equations. Writing ( |l3"| ) explicitely we have, up 
to the coefficients, 

/ /(A,Ai,i/ + l) f(\,» + l,v) /(A-l,//,*/)\ / g(\ + l,n,v) \ 

f(A,fi,v + l) f{\+l,n,v) /(A,/*-l,i/j 9 {\vl + 1,v) 

f(X,fx + l,u) f(\+l,n,v) /(A,m,z/-1) g{\w+l) 

V/(A-l,/i,z/) /(A,/x-1,j/) f{\,p,v-l) / \g{X + X,n+l,v+l)J 

= 0. (15) 

Since this is a homogeneous linear equation the determinant of the coefficient matrix must vanish. For 
simplicity we assume that the matrix is anti-symmetric, that is, the coefficients satisfy ay = — a^j. In 
this case the determinant is the square of Pfaffian, hence the solvability condition turns out to be 

ai4a23/(A + 1,/x, v)f{\ - 1, [i, v) - ai 3 a 2 4/(A, fi + 1, 2/)/(A, /i, - 1, v) 
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+ai2a 34 /(A, n, v + 1)/(A, /x, v - 1) = 0. 



(16) 



In this way the HBDE is reproduced. 

4 Generalization of HBDE to Higher Dimensions 



Our symmetric LBT (|13|) possesses the symmetries of the three dimensional fee lattice space. The fourth 
derivative 84, however, plays a different role in ( |l3|) from others. This owes to the fact that there are 
only three independent variables needed to characterize the four surfaces Si ~ 5*4. This also prevents 
us to extend the space to higher dimensions. 

We like to show, in this section, that the symmetric LBT (|l3| ) can be rewritten in a manifestly 
symmetric form in the four dimensional lattice space if we introduce a parameter s which specifies the 
order of Backlund transformations. Moreover it also enables us to generalize the linear equations to 
arbitrary dimensional lattice space. 

Instead of considering two different kinds of fields we associate the parameter s + 2 to / and s + 1 to 
g and define a new function /(A, fi, v; s) by 

f(\,H,v;s + 2) = f(X,fj,,u), f(X,fi,u;s + l)=g(X ) (i,v). (17) 

We also introduce new variables fcj's with j = 1, 2, 3, 4 by 

A — a — v s u — v — A s v — A — n s 

fei := 1 — , ko := 1 — , kr, := 1 — . 

1 2 -r 4 , 2 2 t 4 , 3 2 T 4 , 

fc 4 := 1 + -. (18) 

It is not difficult to convince ourselves that the operations of dj is equivalent to the operation of dj — d s , 
where dj means d/dkj for all j. Hence the symmetric LBT ( [T^ ) can be written as 

T£ aij (eZL^'-^f) ( e ^f) =0, i = 1,2,3,4, (19) 
i=i 

We notice that ( |l9| ) is totally symmetric for all variables kj, j = 1, 2, 3, 4. 
In terms of the new variables the HBDE turns out to be 

a/(fci + 1, k 2 , k 3 , fc 4 + l)f(ki,k 2 + 1, fc 3 + 1, fc 4 ) 

+Pf(k 1 ,k 2 + l,k 3 ,k 4 + l)f(h + l,k 2 ,k 3 + l,fc 4 ) 

+7/(Ai, k 2 , k 3 + 1, fe 4 + l)/(fci + 1, k 2 + 1, k 3 , k 4 ) = 0. (20) 

The introduction of the fourth variable causes nothing to the contents of the equations. It, however, 
enables us to generalize equations to higher dimensional lattice space. We like to show in the rest of 



this paper that the generalized equations also provide Backlund transformation scheme which generates 
solutions to new nonlinear difference equations. 



For the convenience we define 



Si = e S, 



di - dj 



1=1 



Gtj (/) = ay exp 
Then the n dimensional LBT and its dual are 

n 

J2 G MSj = o, 

n 

J2Ga(S)Sj = o, 



/> G l3 (/) = Oij exp 



Q l + 9 i + 9 3 



L 1=1 



i = 1, 2, 3, • • • , n 
i = 1, 2, 3, • • • , n. 



(21) 
(22) 

(23) 
(24) 



They are the same equation but written differently. The Backlund transformation of the system (|23|) and 
( p4] ) proceeds in a way similar to the LBT of HBDE, which was discussed in §2. For the wave function 
fi to be a solution of the homogeneous equation (E3T) , the potential / must satisfy 



det [G(/)] = 0. 



(25) 



This is a nonlinear difference equation of /. At the same time the solvability condition of (E4) requires 



det G(f) =0 



(26) 



fl25|) and (|26|) are the nonlinear equations to be solved and generalize (|4]). Since ( |25D and ( piq) are not 
the same equation in general, the transformation may not be auto Backlund transformation. In the case 
of n = 4, these two equations are identical and coincide with the square of HBDE, as we have seen above. 
Hence (^5|) and (p6| ) are generalization of HBDE to n dimensional lattice. Notice that at every step of 
the Backlund transformation the value of the parameter s is changed by one. Since s specifies a hyper 
plane in the n dimensional lattice space, the Backlund transformation propagets fields on one plane to 
the next. 

Before closing this paper let us remark briefly the connection of the present results with other works. 
A generalization of HBDE to higher dimmensional lattice space has been discussed in |t3| . In their paper 
the equation itself is bilinear and is satisfied by the same r functions of HBDE. It is a natural extention 
because the space of solutions to the KP-hierarchy itself is symmetric under the choice of three variables 
out of infinite number of variables fci, fe, fe, • • •. We have checked that all soliton type of solutions to the 
higher dimensional HBDE also satisfy our symmetric LBT. Moreover we have proved that there exists 
a solution which satisfies our trilinear equation but does not satisfy any of bilinear Hirota equations 
extended to the four dimensional lattice. 
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An integrablc trilincar partial differential equation (PDE) was discussed in 1 14 , which has a connection 
with the Broer-Kaup system. Generalization to multilinear PDE was also studied in and the Painleve 
properties were examined. The authors claim that no bilinear form corresponding to them has been found. 
Since the equations considered in these references are not differnce but differential equations, we need 
further investigation to see the connection of our work with theirs. 
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